Experimental characterization of quantum dynamics through many-body interactions 
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We report on the implementation of a quantum process tomography (QPT) technique known as 
direct characterization of quantum dynamics (DCQD) applied on coherent and incoherent single- 
qubit processes in a system of trapped ^°Ca^ ions. Using quantum correlations with an ancilla 
qubit, DCQD reduces exponentially the number of experimental configurations required for standard 
QPT. With this technique, the system's relaxation times Ti and T2 were measured with a single 
experimental configuration. We further show the first complete characterization of single-qubit 
processes using a single generalized measurement realized through multi-body correlations with 
three ancilla qubits. 

PACS numbers: 03.67.Ac; 03.65.Wj 



Characterization of quantum dynamics is an impor- 
tant primitive in quantum physics, chemistry, and quan- 
tum information science for determining unknown envi- 
ronmental interactions, estimating Hamiltonian parame- 
ters, and verifying the performance of engineered quan- 
tum devices. This has led to a major effort in develop- 
ing tools for the full characterization of quantum pro- 
cesses, known as quantum process tomography (QPT). 
The standard approach for QPT is resource intensive, 
requiring 12^ experimental configurations for a system 
of N qubits , where each experimental configuration 
consists of the preparation of input probe states and the 
measurement of process outputs. Using ancilla qubits 
but only joint separable measurements, the number of 
experimental configurations is still 12^ [3 5 . However, 
the use of many-body interactions to ancilla qubits in 
the preparation and/or measurements can significantly 
decrease this number to anywhere from 4^ to a single 
configuration depending on the nature and complexity 
of quantum correlations [5 . Using two-body correlations 
DCQD requires up to 4^ experimental configurations for 
full quantum process tomography, and in particular only 
one experimental setting for estimating certain parame- 
ters (e.g. relaxation times) [6l[7] (see Ref. [3 [9] for partial 
and non-scalable implementations of DCQD). 

Alternative tomographic methods such as randomized 
benchmarking [10 , selective and efficient QPT pT| [T2] 
and compressed sensing for quantum process tomogra- 
phy [TSHIS] have recently been developed to overcome 
the exponential increase of the required experimental 
configurations. Generally, these methods are tailored to 
estimate a polynomial number of effective parameters, 
such as gate fidelity [10 or when we can make a sparse 
quantum process/Hamiltonian assumption from a priori 
knowledge about the quantum system [15]. For exam- 
ple, the estimation of the dynamical parameters Ti and 



T2 (longitudinal and transverse relaxation times |T) is a 
task involving two non-commuting observables (e.g. cFx 
and (Jz) that cannot be measured simultaneously. These 
parameters describe the influence of noise on atomic-, 
molecular- and spin-based systems induced by the inter- 
action with the environment. An alternative approach 
based on DCQD, henceforth called Direct Characteriza- 
tion of Relaxation Times (DCRT), enables the measure- 
ment of both Ti and T2 simultaneously with a single 
experimental configuration [16 . 

In this work, we apply the DCQD technique and 
extensions on a system of trapped ^^Ca+ ions. Single- 
qubit processes are reconstructed with four experimental 
configurations using DCQD, and alternatively with just 
a single configuration using a generalized measurement 
(GM). In addition, we quantify the relaxation times 
Ti and T2 in our system with a single configuration. 
This technique can also characterize more realistic 
environments affecting not only the probe but also the 
ancilla qubit collectively. 

Every process acting on a quantum system can be de- 
scribed by a complete positive map E mapping the input 
state p onto the output state p' . For a single qubit this 
can be written as 

4 

E \ p^ p' = ^ Xm,n CTm P Cr^, (1) 

m,n—l 

with cr^, an the Pauli operators {1, a^^, a^, a^} and x 
a semi-positive matrix containing complete information 
about the process. In standard quantum process tomog- 
raphy (SQPT) the process is applied to four input states 
and followed by full state tomography of each output 
state, which for a trace preserving map consists of three 
measurements, resulting in 4 x 3 = 12 experimental con- 
figurations. In DCQD these four input states are replaced 
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TABLE I. Input states and BSM basis used for the imple- 
mentation of DCQD = i^^, = ^^^P)- The 
determinant of A in Eq. (2) is maximized for a = cos(^ 
and /3 = exp(i|^) sin(^)) to ensure the invert ibility ^17^ . 
The BSM is realized by a measurement with the projectors 



P. = {|$±)($±|,|vl/±)(vl/±|}. 



FIG. 1. (color online) Procedure to characterize a single-qubit 
process with DCQD and a GM. In DCQD (a) each exper- 
imental configuration consists of the preparation of one of 
four input states pj entangled between the system ion S and 
the ancilla ion A. The process E is applied on S followed by a 
BSM on the output state E{pj), which consists of a single MS 
operation followed by a projection onto the computational ba- 
sis, (b) Generalized measurement via many body interactions 
(see text). 

by four entangled states between the system qubit S and 
an ancilla qubit A, and the state tomography is replaced 
by a single Bell-state measurement (BSM), as shown in 
Fig. [ija), with a total of 4 x 1 = 4 experimental con- 
figurations (Bell states |^^) and |<l>^) as defined in Ta- 
ble [l|. The probabilities pij of measuring the Bell-state 
projector Pi for each input state pj shown in Table [l| are 
determined, according to Refs. [6l[T7], by 

4 

= TriPi E{pj)) = ^ Xm,n A^,n (2) 
m,n=l 

Therefore the process matrix x can be calculated directly 
by linear inversion of the matrix A. The set of input 
states pj and the Bell-state projectors Pi have to be de- 
termined such that the 16 equations in Eq. (2) are linearly 
independent, which ensures that A is invertible (Table[l|. 

Our experiments were realized on a system consist- 
ing of ^^Ca+ ions confined to a string in a linear Paul 
trap [18 . Each ion represents a logical qubit which is en- 
coded in the electronic levels D^/2{^ = ~l/2) = |0) and 
S'i/2(m = —1/2) = |1). Each experimental cycle consists 
of an initialization of the ions in their internal electronic 
and motional ground states followed by a coherent ma- 
nipulation of the qubits and finally a detection of the 
quantum state. State initialization is realized by optical 
pumping into the 6'i/2(m = — 1/2) state after cooling the 
axial centre-of-mass mode to the motional ground state. 
The manipulation of the qubits is implemented by co- 
herently exciting the S1/2 ^ ^5/2 quadrupole transition 
with laser pulses. Finally, the population of the qubit 
states is measured by exciting the 51/2 ^ P1/2 transi- 



tion and detecting the fluorescence light, using electron 
shelving [19]. Our setup is capable of realizing collec- 
tive qubit rotations U{6^(j)) = exp(— z| ^Jsin(^)(j^*^ + 

cos{(j))ax^]) via a laser beam addressing the entire register 
as well as M0lmer-S0renson entangling gates MS{0^ (j)) = 
exp(-if [EiSin(0)4'^ +cos(^)cri'^]2) [20l[21]. Addition- 
ally we are able to perform single-qubit rotations on the 
i-ih ion of the form U^\0) = exp(— z|cri*^) by an off- 
resonant laser beam, which addresses individual ions. 
The input states for DCQD of Table I are prepared by 
applying collective entangling operations and qubit ro- 
tations as shown in Fig. [ija). For example, the input 
state p2 is created by the non-maximally entangling op- 
eration M5'(J,7r). Our two-qubit entangling operation 
generates Bell states with a fidelity of ^ 99% in 120 ps. 

The BSM is experimentally realized by a max- 
imally entangling operation MS{^^ ^), which maps 
from the Bell-state basis to the computational ba- 
sis {|00), |01), |10), |11)}, followed by individual-ion- 
resolving fluorescence detection with a CCD camera. For 
example, consider the first input state pi = |<l>+)(<l>+|. If 
the process E is the identity 1, the expectation value of 
the BSM-projector Pi is 1 which is equivalent to detect- 
ing both ions in the state |11) after the BSM. If a bit flip 
occurs on the system ion, the output state is then mapped 
onto the state |01) by the BSM ((P2) = 1). The consid- 
erations are similar for a phase-flip, or bit- and phase-flip 
processes. Therefore, the diagonal elements Xm,m of the 
superoperator x corresponding to 1, (Jy^ and are 
detected by a single input state in combination with one 
BSM. 

We demonstrate the DCQD method by characterizing 
the full quantum process of implemented unitary rota- 
tions and cFy as well as non-unitary processes such 
as amplitude- and phase-damping ^ . The x matri- 
ces reconstructed from the measured probabilities, are 
shown in Fig. |2|a,b) for (Jx and CFy rotations. A single- 
qubit process can also be visualized by transforming the 
pure states lying on a Bloch sphere. In this Bloch sphere 
representation, decohering processes map the unit Bloch 
sphere (shown as a transparent mesh) to an ellipsoid of 
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smaller volume [T]. Implemented amplitude- and phase- 
damping processes taking place with a 60% probability 
are shown in this representation in Fig. [2|d,f) [22]. For 
each input state the experiment was repeated up to 250 
times for statistical averaging. All processes were recon- 
structed with a maximum likelihood algorithm to ensure 
trace preservation and positivity of the superoperator 
X ^3 . The fidelity F of a process describes the overlap 
between the measured Xmeas and the ideal superoperator 
Xid' For each process we calculate the overlap between 
Xmeas and Xid usiug the Uhlmann-Jozsa fidelity with the 
Jamiolkowski isomorphism [23l |24] . Table |ll| shows the 
calculated fidelities for the implemented DCQD and for 
SQPT. The uncertainty in the fidelity was estimated by 
parametric bootstrapping based on projection noise in 
our measurement [25] . 



Target process 


DCQD, F (%) 


SQPT 


, F (%) 


1 


97.5 ±0.6 


98.1 


± 1.3 


CTx 


96.5 ± 1.0 


98.1 


± 1.3 


ay 


96.6 =b 1.4 


97.5 


± 1.4 


amplitude damping 


95.3 ± 1.9 


95.2 


±2.7 


phase damping 


97.4 ±0.8 


95.7 


±0.8 



TABLE II. Calculated process fidelities F between imple- 
mented and target processes as characterized with DCQD and 
SQPT. All processes were measured with a total of 1000 ex- 
perimental cycles, which correspond to 1000/4 cycles per ex- 
perimental configuration for DCQD and 1000/(4x3)^84 for 
SQPT. The SQPT of the phase damping process was mea- 
sured with a total of 3000 experimental cycles. 



Full QPT of a single-qubit process is also possible with 
a single experimental configuration by using additional 
ancillas and a generalized measurement (GM). Here, we 
expand the dimension of the Hilbert space Ha^Hs, with 
the system Hilbert space Hs and the ancilla Hilbert space 
Ha^ such that the dimension of the total Hilbert space is 
equal to the number of free parameters in the process ma- 
trix X [5^ . For a single-qubit process one has to determine 
all 16 superoperator elements Xm,n which leads to an 8 di- 
mensional ancilla Hilbert space. Therefore we used three 
ancilla qubits Ai, A2 and As to quantify a full process E 
acting on the system qubit S. This GM is realized by en- 
tangling the system and ancilla qubits using many-body 
interactions [20^, '2T , then applying the process £^ on S 
and finally performing BSMs on two pairs. Figure [ijb) 
shows the sequence implemented for this GM which pro- 
ceeds as follows. First, we create an entangled input state 
using maximally and non-maximally entangling M0lmer - 
S0renson interactions in combination with global and ad- 
dressed single-qubit rotations. After applying the process 
£^ on S we perform a pairwise BSM on the combined out- 
put state by implementing two non-maximally entangling 
operations MS{j) and two addressed AC-Stark pulses 

U^\7r) and [/^^^(tt), which separate the entangled sys- 




FIG. 2. (color online) Experimental results of DCQD for uni- 
tary and decoherence processes, (a-b) Results of the mea- 
sured superoperator x foi" the rotation operations [/(tt, 0) in 
(a) and [/(7r,7r/2) in (b). Ideally, the target processes have 
only nonzero elements at positions indicated by the orange- 
bordered bars, (c-d) Bloch sphere representation of the ideal 
(c) and measured (d) amplitude damping process with 60% 
probability [22]. (e-f) Bloch sphere representation of the ideal 
(e) and measured (f) phase damping process with 60% proba- 
bility [22] . Bloch sphere axes in black evolve into the spheroid 
primed axes in blue. A slight imperfection due to residual 
light on the ancilla ion can be observed as a rotation of the 
spheroids in the measured decohering processes. 



tem H{S^ Ai, A2, A3) into a product state of two subsys- 
tems H{Ai, As)^ H{S^A2). These operations are equiv- 
alent to two pairwise maximally entangling gates MS{^) 
acting on the two subsystems H{Ai^As) and i^(5', A2). 
The 16 results of the measurement are directly linked to 
the 16 superoperator elements Xm,n by a matrix A sim- 
ilar to Eq. (2). Using this technique we reconstructed 
unitary processes {t^dx = [/(7r,0),cr^ = [/(tt, |),cr2 = 

U^z^\7r)} acting on a single qubit with a fidelity of 
{99.70 ± 0.02, 97.30 ± 0.29, 99.80 ± 0.01, 99.40 ± 0.02}%. 
All processes were measured with a total of 5000 cycles. 

In contrast to previous QPT measurements of engi- 
neered processes, the process of phase (amplitude) damp- 
ing occurs naturally in our system due to magnetic field 
fiuctuations (spontaneous decay) [26]. The dynamical 
parameters Ti and T2 can, however, be determined si- 
multaneously with only the first input state pi being sub- 
ject to the DCQD scheme even if the damping processes 
act collectively on both qubits (as in our experimental 
system [26 ). This method, named DCRT above, consists 
of preparing an input Bell state pi = |<I>+)($+|, exposing 
both qubits to the damping processes for a time t and a 
final BSM, which yields the diagonal elements Xi,i of the 
process matrix. As described in the supplementary ma- 
terial [22^ and assuming Markovian noise the dynamical 
parameters are then given by 
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e "^2* = xi^i - X4,4 (3) 
= Tr{[|^+)($+|-|$-)($-|]E(|$+)($+|)}, 

_ 2t_ t_ 

1 + 2e ^1 - 2e ^1 = 1 - 2(x2,2 + X3,3), (4) 

with N the number of ions. From the entries of the x 
matrix corresponding to 1 and {(Jx and (Ty) opera- 
tions, T2 (Ti) depends on the probabiUty that no error 
or phase flips (bit flips) occur on the entire system. A 
flt of DCRT measurements Xi,i to Eqns. (3-4) at different 
times t thus yields Ti and T2 using a single experimental 
conflguration. We explored this DCRT technique in our 
experimental system. The measurement results of the de- 
coherence estimation are shown in Fig. |3ja). The green 
dots show the difference between the diagonal elements 
Xi,i and X4,4 as a function of the waiting time t. The 
spontaneous decay of the system is shown in Fig. [sj^b) 
by plotting 1 — 2(x2,2 + X3,3) as a function of time. For 
every waiting time t the experiment was repeated up to 
250 times to gain signiflcant statistics. 

We can compare the DCRT technique with two 
traditional methods that use product input states: 
Ramsey-contrast measurements for phase-decoherence 
estimation and direct spontaneous-decay measure- 
ments [27]. A Ramsey- contrast measurement is realized 
by initializing the ion in the state (|0) + |1))/V2 by 
a global rotation [/(|,0), followed by a waiting time 
t and flnally applying a second rotation [/(|,^) in 
which the phase (j) is varied. The observed contrast 
as a function of (j) corresponds to the preserved phase 
coherence. Spontaneous-decay measurements, instead, 
consist of measuring the probability of detecting both 
ions in the excited state |0) as a function of time. The 
results of these Ramsey-contrast (spontaneous-decay) 
measurements are shown in Fig. [sja) (Fig. jsjb)) as red 
diamonds (blue triangles). The measured relaxation 
times corresponding to the traditional methods are 
called Tf^^ and T^^"""^ . The exponential flt (green 
line) of Eq. (3) to the data was estimated with N = 2 
(collective dephasing) and yields T2*^^-^=18.8(5) ms. 
The Ramsey-contrast measurements (red diamonds) 
were carried out on a single ion and yield a coherence 
time of T|^^^=19.4(8) ms. The green dotted line in 
Fig. |3ja) corresponds to the fltted function to Eq. (3) 
(green line) with = 1 instead of A" = 2 and shows 
good agreement with the single-ion Ramsey-contrast 
measurement. Therefore the DCRT technique enables 
the characterization of the phase decoherence of the 
collective system (green line) and also gives a conclusion 
about the phase decoherence of a single ion (green 
dotted line). An exponential flt of the decay data of 
Fig. Isjb) to Eq. (4) gives the characteristic lifetime 
rj^Dckr ^ 1130(47) ms for the DCRT technique (green 



line) and Tf = 1160(30) ms for the traditional method 
(blue dotted line), which are in good agreement with 
previously measured values [28 of 1148(18) ms. 

In summary, we have experimentally demonstrated two 
different approaches for the full characterization of single- 
qubit quantum processes, lowering the required experi- 
mental conflgurations from 12 to 4 using DCQD and a 
single conflguration via the CM method. The reconstruc- 
tion of coherent and incoherent processes was shown with 
fldelities of > 97% using DCQD. In particular, we have 
observed a lower statistical uncertainty of the fldelity of 
some of the processes compared to the SQPT. Never- 
theless, a matter of further investigations is a compari- 
son of the scaling in the number of experimental cycles 
required for the SQPT and DCQD to achieve a target 
uncertainty in the fldelity (e.g. see identity process in 
Table |ll|. The DCRT technique, based on the DCQD 
protocol, was used as a powerful tool to characterize the 
noise in our system by measuring the relaxation times 
Ti and T2 simultaneously with one experimental set- 
ting. This technique indicates good agreement with tra- 
ditional methods as Ramsey-contrast and spontaneous 
decay measurement. In principle, there is an improve- 
ment of a factor of two in the measurement time if Ti is 
of the same order of magnitude as T2, which is not the 
case for our setup. In contrast, spin-based solid state sys- 
tems are collectively affected by noise and Ti^T2, which 
would lead to a signiflcant improvement of the measure- 
ment time [29 . Another application of DCRT could be 
for biological systems where dissipative dynamics play a 
crucial role [30l [31] . The same measurement procedure 
can also be used as a tool to quantify Hamiltonian pa- 
rameters efliciently [71 116] . 
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Appendix 



Implementation of amplitude- and phase-damping 



For the complete quantum process tomography using 
DCQD it was considered that the process of interest only 
acts on the system qubit. In general due to residual 
light on the ancilla ion (crosstalk), it is experimentally 
challenging to apply single-qubit operations on one ion 
without affecting the second one. The unitary opera- 
tions on the system ion S without affecting the ancilla ion 
A is experimentally implemented by the following refo- 
cusing technique: (i) Applying a global rotation [/(|,^) 
around the X or Y" axis of the Bloch sphere, (ii) Shifting 
the phase of the second ion by tt via an AC-Stark pulse 

(2) 

(tt). (iii) Repetition of the first step. This effectively 
reduces the crosstalk, for single qubit operations roughly 
by a factor of 6, ideally to less than 1%. 
The process of phase damping is implemented by the 
following steps (see Fig.[4][a)): (i) Hiding the population 
of the 1^5/2 (^ = —1/2) state of both ions by applying a 
TT-pulse on the 1)5/2 (^ = —1/2) ^ Si/2{m = 1/2) tran- 
sition, (ii) Transfering a certain amount of the popula- 
tion of the Si/2{^ = —1/2) state to the D^i2{m = —5/2) 
state. The excitation probability is experimentally con- 
trolled by the pulse length t and is given by sin(7)^, with 
7 = ^ and T the Rabi-oscillation period time, (iii) Re- 
pumping the transfered population into the P3/2 state 
by a laser pulse at a wavelength of 854 nm whereupon 
the system qubit decays spontaneously into the ground 
state, (iv) Finally the second step is repeated to reverse 
the hiding process. Because of the spontaneous decay 
after the third step the state loses its phase information. 
The hiding process described in the first step is necessary 
to prevent the Dr^i2{m = —1/2) state from not being af- 
fected by the repumping process which would have the 
same effect as amplitude damping. 

Amplitude damping is carried out within two steps 
(see Fig. |4]^b)). First a certain amount of population is 
coherently transfered from the D^i2{m = —1/2) to the 
Sii2{Tn = 1/2) state by an addressed laser pulse analog to 
step (ii) of the phase damping process. Then the whole 
population of the Si/2{m = 1/2) state is transfered to 
the Si/2{m = —1/2) ground state by an optical pump- 
ing process performed by a cr~ -polarized laser beam at 
397 nm exciting the ion first to the Pi/2(^ = —1/2) state 
followed by spontaneous decay to the Si/2{m = —1/2). 
The ancilla ion is not affected by the optical pumping 
process because the cr~-polarized pulse can only excite 
the Sii2{Tn = 1/2) ground state. 



Direct characterization of collective longitudinal and 
transverse relaxation processes 

Let us consider a global quantum homogenization or 
thermalization process acting on some qubits for a time t. 
Thus for each single-qubit density matrix p, with poo = ct 
and poi = ^ in the computational basis, evolves into the 
state p{t) with poo{t) = {a — ao) exp(— t/Ti) + ao and 
Poi{t) = 6exp(— t/T2), where Ti and T2 (T2 ^ 2Ti) are 
relaxations and dephasing time-scales of the system, re- 
spectively. That is, the system approaches an equilibrium 
state identified by ao G [0, 1]. 

In order to measure these time scales using the DCQD ap- 
proach we need to create quantum correlations between 
each pair of neighboring qubits, e.g., in the form of a 
Bell-state p = |^+) (^+| where |^+) = (|00) + |11))/V2, 
and then let all of the qubits evolve for a time t lead- 
ing to state S{p). Thus by performing a Bell- state mea- 
surement between the same neighboring qubits we obtain 
(ESI, EH): 



Tr[P^^S{p)] - Tr[P^^£{p)] = e"^ 



where P^^' = \B^) {b^' , and \B^) for k = 1,2,3,4 cor- 
responds to the Bell-states |^+), |^+), |^"), and |$"), 
respectively, where = (|00) ± |11))/a/2, = 

(|01) ± |10))/V^. Using the fact that Tr[P^^£^{p)\ = 
Xi,i and Tr[P^^8^ {p)] = X4,4 we obtain: 



Xl,l - X4,4. 



(5) 



It is important to mention that this expression is only 
valid for Markovian noise acting on the entire system. For 
our system it was observed [26] that the phase decoher- 
ence of Greenberg-Horne-Zeilinger (GHZ) states scales 
with exp(— A/'^t) instead of exp(— TVt), with N the num- 
ber of ions. Therefore Eq. (1) has to be modified, which 
leads to the following result for the phase-decoherence: 



Xl,l - X4,4. 



(6) 



In the case of Ti, the outcome of a BSM within the 
DCQD formalism yields: 



1 - 2{Tr[P''8{p)] + Tr[P-''£{p)]) = 
= (1 - 2aof{l - 2e-'T) + (2 + 4ao(ao - l))e"t . (7) 

For our system the equilibrium state is described by the 
ground state and therefore ao = 1, which leads to: 



1 - 2(X22 + X33) = (1 - 2e ^i) + 2e 
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FIG. 4. (color online) Schematic of the phase damping (a) and amplitude damping (b) processes acting on the system ion S: 
(a) Transfer of the Si/2(rn = —1/2) population with a probability of sin(7)^ to the P3/2 state followed by optical repumping to 
the ground state while the population of the D^/2{m = —1/2) state is hidden in the 5*1/2 (m =1/2) state, (b) The amplitude 
damping process is carried out by coherently transfering the population of the 1)5/2(771 = ~l/2) state to the Si/2{'rn =1/2) 
state of the system ion followed by optical pumping to the ground state Si/2(rn — —1/2). The ancilla ion is not affected during 
the whole processes, as shown in the inset figures. 



For a unital quantum homogenization process {E{I) = /), 
we have ao = ^ (i.e., a completely stochastic equilibrium 
state). Thus the relation for Ti becomes: 

- — 

e ^1=1- 2(X2,2 +X3,3). 



It is remarkable that even if we use the relations devel- 
oped for the ideal DCQD scheme, assuming Ti and T2 



acting only on the system qubit (considering the other 
system as a reference signal or a noiseless ancilla), we 
can still obtain both T2 and Ti only to be smaller than 
the actual value by a factor of 2. Note that, due to or- 
thogonality of the BSM outcomes, it is easy to unambigu- 
ously distinguish Ti from T2. Traditionally, in order to 
measure the longitudinal and transverse relaxation times, 
one needs to measure two non- commutative observables 
(e.g., Pauli operators and ax) on two sub-ensembles 
of identical systems. 



